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SUMMARY 


Several  basic  conditions  which  determine  the  shape  of  the  free  fluid 
surface  due  to  the  vertical  immersion  of  an  infinite  wedge  in  a  perfect 
fluid  are  analyzed.  Continuity  of  the  fluid  flow,  dynamic  similarity,  and 
finally  the  conditions  for  irrotational  potential  flow  are  applied  to  es¬ 
timated  surface  shapes  until,  by  an  iterative  solution  method,  a  unique 
surface  shape  is  obtained  for  a  given  wedge  angle.  This  method  of  analysis 
is  carried  out  for  several  wedges  (deadrise  angles  up  to  50°),  and  the 
corresponding  free  surface  shapes  are  constructed  in  detail.  Empirical 
formulae  for  the  relative  spray  thickness  are  also  given. 

Knowledge  of  the  free  surface  shape  enables  a  solution  for  the 
velocity  and  potential  everywhere  on  the  free  surface,  and  also  in  the 
wedge  itself.  From  these  data,  the  impact  force  and  pressure  distribution 
on  the  wedge  are  computed.  Comparison  of  the  pressures  with  a  previous 
simpler  analysis  based  on  an  analogy  with  the  flow  around  an  expanding 
flat  plate  shows  reasonable-  agreement  at  low  angles,  but  a  large  percent¬ 
age  reduction  in  pressure  at  high  deadrise  angles. 

The  study  wasmade  at  the  Experimental  Towing  Tank,  Stevens  Institute 
of  Technology,  under  the  sponsorship  of  the  Office  of  Naval  Research, 
Department  of  the  Navy. 


INTRODUCTION 


This  report  is  one  of  a  series  prepared  at  the  Experimental  Towing 
Tank,  Stevens  Institute  of  Technology,  in  connection  with  the  research  on 
planing  surfaces  conducted  under  Contract  No.  N6onr-247,  Task  Order  IV, 
with  the  Office  of  Naval  Research.  The  research  project,  which  includes 
theoretical  and  experimental  investigations  on  the  problem  of  planing  on 
the  surface  of  water,  is  divided  into  two  general  groupings: 

a)  Investigation  of  the  fundamental  nature  of  the  hydrodynamic 
planing  process. 

b)  Collection  and  organization  of  seaplane  and  flying-boat  design 
data  for  the  establishment  of  rational  design  methods  and  for 
comparison  with  theoretical  results. 

Group  (a)  is  primarily  of  a  theoretical  or  mathematical  nature, 
dealing  with  the  detailed  fluid  flow  in  the  vicinity  of  a  planing  surface. 
Results  of  these  theoretical  studies  are  presented  in  Refeiences  1,  2, 
and  3,  with  additional  reports  in  the  process  of  preparation.  Group  (b) 
is  primarily  experimental  and  empirical,  dealing  with  an  over-all  evalu¬ 
ation  of  the  effects  of  the  many  parameters  that  influence  the  planing 
process.  Results  of  the  empirical  studies  completed  to  date  are  presented 
in  References  4,  5,  6,  and  7. 


The  present  report  is  a  portion  of  the  investigation  into  the  funda¬ 
mental  nature  of  the  hydrodynamic  planing  process.  It  deals  in  detail 
with  t..e  rise  of  the  surface,  formation  of  spray,  pressures,  and  loads 
'"hich  occur  during  the  normal  penetration  of  a  fluid  surface  by  an  in¬ 
finite  wedge. 

The  general  appearance  of  the  spray  and  wave  surface  which  results 
from  the  penetration  of  a  fluid  surface  by  a  wedge  shape  has  been  experi¬ 
mentally  noted  by  many  observers.  However,  the  transitory  nature  of  the 
problem  and  the  small  measurements  involved  have  made  accurate  experimen¬ 
tal  data  difficult  to  obtain.  Loads  and  pressure  distributions  for  wedge 
immersions  have  been  determined  experimentally  (because  of  the  importance 
of  the  problem  in  seaplane  design),  but  only  at  the  expense  of  intricate 
and  costly  test  procedures. 

Analysis  and  understanding  of  the  experimental  data  obtained  have 
be;n  hampered  by  •>  lack  of  a  basic  analytical  solution  to  the  penetration 
problem.  The  earlier  work  of  Wagner  (Reference  8),  brief  and  cryptic  as  it 
is,  has  nevertheless  provided  a  groundwork  for  the  theoretical  study  of 
the  many  problems  associated  with  planing  and  impact.  Previous  reports  of 
this  series  have  amplified  and  extended  Wagner’s  analyses  for  the  detailed 
flow  in  the  spray-root  and  trail ing-edge  regions  and,  by  the  expanding- 
plate  analogy,  for  the  loads  and  pressures  on  a  wedge  penetrating  the 
fluid  surface.  These  analyses  were  limited  to  relatively  low  angles  by 
the  assumptions  required  for  their  solutions.  The  removal  of  the  restric¬ 
tion  to  low  deadrise  angle  would  make  the  results  of  these  theoretical 
studies  of  greater  practical  engineering  value. 

It  is  the  purpose  of  this  report  to  present  the  more  general  method 
of  analysis  suggested  by  Wagner  which  yields  the  free  surface  shape  and 
associated  pressures  and  loads  during  the  symmetrical  penetration  of  a 
fluid  surface  by  a  wedge  of  any  angle. 


SYMBOLS 


a-b  points  on  free  surface  of  fluid;  also  distance  of  points  from  origin 
o 

•'  nominal  wetted  half-width 

d  subscript  pertaining  to  doublet 

id  subscript  pertaining  to  image  doublet 

k  U/V  at  point  on  surface 

n  distance  normal  to  free  surface,  *  inward 

0  origin  of  x,y  axes 

o’  edge  of  spray,  origin  for  s 

P  pressure  at  any  point  in  flow  field  [-pCdi/dt)  -  p(U2/2)  +  pF(t)) 

R  distance  to  doublet  center 

r  distance  from  origin  to  point  in  fluid  (usually  on  surface  curve) 
s  distance  along  surface  from  its  origin  at  solid  boundary  (Figure  2); 
also  general  point  on  free  surface  or  wedge 


f  time  _ 

U  resultant  velocity  at  any  point/  n2  +  v2 

u  velocity  component  ||  x  axis;  1  ~  ;  d<fi/dx,  '&fr/dy 

V  velocity  of  wedge  or  field  velocity  at 

v  velocity  component  II  y  axis;  +  f  >  <W3y#-<W3x 
x  horizontal  distance  to  <£.  symmetry 

x  position  of  peak  pressure  on  bottom 

y0  vertical  distance  to  apex  of  wedge  (from  undisturbed  water  surface) 

a  angle  betv/een  velocity  vector,  U,  and  radius,  r 

0  angle  of  deadrise  (as  a  subscript,  denotes  polar  coordinate  system 
(Figure  ,)) 

y  angular  orientation  with  respect  to  doublet 

?  spray  thickness 

4  angle  between  the  radial  velocity  and  the  surface  velocity  past  the 
radius 

V  y  coordinate  of  particular  fluid  particle  (especially  surface  par¬ 

ticles) 

&  angle  between  radius  from  doublet  and  free  surface  at  any  point 

p  angle  subtended  at  the  doublet  by  the  line  yQo‘ 

x  coordinate  of  particular  fluid  particle  (especially  surface  par¬ 
ticles) 

p  mass  density  of  fluid 

0  velocity  potential 

0  stream  function 

V  differential  operator  (V0  =  grad0) 

V2  Laplacian  operator  <=  div  grad 


THEORETICAL  ANALYSIS 


CONDITIONS  OF  THE  PROBLEM 

The  problem  is  treated  as  a  two-dimensional  one,  and  the  infinite 
wedge  is  assumed  to  be  symmetrical  about  the  vertical  centerline  plane, 
and  perfectly  rigid.  For  this  first  analysis,  the  velocity  of  penetration 
is  taken  as  a  constant  (very  high  mass  or  forced  penetration). 

The  fluid  is  assumed  incompressible  and  frictionless.  The  flow  field 
is  assumed  to  be  initially  irrotational ,  and  thus  must  remain  so  under 
the  first  assumption  of  the  perfect  fluid.  The  velocity  of  penetration  is 
assumed  to  be  sufficiently  high  and  the  size  of  the  wedge  great  enough 
that  the  effects  of  gravity  and  surface  tension  may  be  safely  neglected. 

The  initial  penetration  of  the  wedge  splits  the  previously  undis¬ 
turbed  fluid  surface  into  two  surface  sheets  On  the  basis  of  the  assump¬ 
tions  of  perfect  fluid  and  irrotational  flow,  it  follows  that  particles 
on  the  original  surface  remain  and  no  new  particles  are  added  (Reference 
9)  and,  since  the  fluid  is  incompressible,  the  surface  as  deflected  by 
the  wedge  penetration  must  have  a  constant  length  of  arc. 

These  conditions  correspond  to  the  requirements  for  a  potential 


flow.  For  convenience  in  application  to  the  problem,  however,  three  con¬ 
ditions  are  abstracted  which  define  the  shape  of  the  free  surface: 

1)  Continuity:  The  fluid  displaced  by  the  wedge  must  appear  above 
the  original  water  surface  in  the  form  of  a  wave  and/or  spray 
(incompressibility),  and  the  arc  length  of  the  surface  must  be 
constant. 

2)  Similarity:  Since  the  wedge  is  symmetrical  and  the  immersion 
is  normal  to  the  undisturbed  surface,  the  field  of  flow  is  di¬ 
vided  into  two  halves  by  the  centerline  of  symmetry.  The  neglect 
of  gravity  forces  and  scale  effects  leads  to  the  definition  of 
surface  shape  and  dynamic  state  of  the  field  of  flow  entirely 
in  terms  of  the  wedge  penetration  (distance  and  velocity).  At 
constant  penetration  velocity,  the  flow  field  and  surface  shape 
must  then  be  geometrically  and  dynamically  similar  at  all  times. 

3)  Irrotationality :  The  fluid  flow  field  may  be  represented  by 
the  potential  function,  <t>.  For  irrotational  flow,  the  relation 
V2d>  =  0  must  be  satisfied  everywhere  in  the  flow  field,  where 
V2  is  the  Laplacian  operator  Jiv  giad. 


METHOD  OF  APPLICATION  OF  GENERAL  CONDITIONS 

CONI INU'TY 

There  are  any  number  of  possible  solutions  for  the  free  surface 
shape  on  the  basis  of  continuity  alone.  Two  extreme  solutions  are  shown 
in  Figure  1:  one  limiting  case  is  all  wave  with  infinitesimally  thick 
spray,  while  the  other  is  the  opposite  --  all  spray  and  no  wave.  Between 
these  extremes  lie  the  combinations  of  various  spray  and  wave  shapes  as 
indicated  in  Figure  2- 

One  analytic  expression  for  the  ‘all  wave’  condition  is  obtained 
from  the  equations  of  flow  about  a  flat  plate  as  given  in  References  3 
and  8.  From  consideration  of  the  potential  flow  about  a  flat  plate,  with 
the  flat-plate  width  treated  as  a  variable,  it  is  possible  to  determine 
the  velocity  distribution  and  then  the  effective  surface  rise  in  the  plane 
of  the  plate  (as  shown  in  Figure  1).  This  solution  gives  an  analytic  ex¬ 
pression  for  the  wave  rise,  7?.  The  rate  of  expansion  of  the  plate  deter¬ 
mines  the  equivalent  deadrise,  to  which  this  analytic  solution  is  as¬ 
sumed  to  apply.  It  is  simple  to  demonstrate  by  integration  of  the  area 
under  the  wave  that  the  continuity  condition  is  satisfied. 

From  Figure  1,  the  areas  A  and  B  for  thi  'all  wave'  condition  are 

>1  .  MM.fl.2j  (1) 

and 
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B  =  _t^i£  (c  _JLc)2+  /n£fx 


(2) 


From  Equation  38  of  Reference  8,  (or  Equations  11  and  13  of  Reference  3), 
the  wave  rise  q  is  given  by: 


'/?  =  itan£  x  arcsin  £  -  c 

77  XV 


(S) 


”  CO 

f  T)  dx  -  c2  /^(-~  arcsin  —  - 

*:  2  y  \vc2  x  vc I 


-  c2  [~^£_  arcsin 

^  L  c? 


£  +  1  A2  .  c2  ~il 

*  nc  Y  '/7cjc 


x-<» 


=  tan/3  ^  2 


R'-i)]  ■ 


(4) 


Therefore,  combining  Equations  (2)  and  (4)  gives: 

B  =  tan^  c2  (l  -  1)  2  +  c2  (A  -  A 

2  \  vl  2  Vtj  / 


tanl  c  2 


(4)*  ■ 


(5) 


which  is  equivalent  to  the  expression  for  the  area  A. 

The  ‘all  spray’  condition  could  be  solved  directly  by  simultaneous 
equations  for  the  area  and  length  equalities  required  by  the  continuity 
condition: 


B 


and  oa  =  o  e 


(6) 


However,  these  two  cases  are  evidently  not  real  possibilities  because  of 
the  sharp  fluid  corner  which  would  be  required  either  at  the  wedge  or  at 
the  still  water  surface.  They  merely  represent  the  extremes  between  which 
an  infinite  variation  of  wave  and  surface  shape  may  occur. 

In  general,  any  thickness  of  spray  may  be  associated  with  an  un¬ 
limited  variety  of  waves  and  still  not  violate  the  continuity  requirement. 

SIMILARITY 

The  principle  of  similarity  for  the  proportional  increase  of  the 
wave  and  spray  with  continued  penetration  is  the  second  basic  considera"^ 
tion.  Not  only  must  the  displaced  fluid  appear  above  the  initially  undis- 


turbed  surface,  but  the  wave  and  spray  must  be  increasing  in  such  a  man¬ 
ner  that  the  shape  and  velocities  are  everywhere  determined  by  the  posi¬ 
tion  and  penetration  velocity  of  the  wedge.  By  a  geometrical  construction 
and  integration  process  proposed  by  Wagner  (Reference  8),  it  is  possible 
to  relate  the  surface  shape  and  velocities  so  that  for  any  spray  thick¬ 
ness,  one  corresponding  wave  and  spray-root  shape  which  meets  the  require¬ 
ments  of  similarity  and  continuity  is  determined. 

In  an  ideal  fluid  where  gravity,  viscosity,  and  surface  tension  ef¬ 
fects  are  neglected,  the  shape  of  the  free  surface  is  uniquely  determined 
by  the  shape  of  the  body,  and  the  size  of  the  shapes  is  determined  by  the 
amount  of  penetration.  Thus,  fora  simple  wedge  penetration,  the  free  sur¬ 
face  siiape  is  always  the  same. 

For  the  same  relative  positions  in  the  fluid  and  on  the  surface, 
the  dynamic  conditions  of  pressure  and  velocity  depend  upon  the  surface 
shape  and  must  bear  the  same  relationship  to  the  penetration  velocity  at 
every  instant  of  time.  It  should  be  noted  that  this  statement  applies  to 
a  relative  position  in  the  fluid  and  not  to  a  particular  fluid  element 
since,  as  it  will  be  shown  later,  particular  particles  move  from  one  posi¬ 
tion  to  another  and  are  subject  to  varying  pressure  and  velocity  with  the 
passage  of  time 

For  example,  consider  the  several  stages  of  the  penetration  process 
illustrated  in  Figure  3,  Initial  contact  is  made  at  point  o,  splitting 
the  surface  layer.  Penetration  into  the  fluid  raises  the  initially  flat 
surface  into  a  wave  and  deflects  the  nearby  particles  out  along  the  wedge 
surface  as  spray.  (The  transition  area  where  the  wave  surface  is  rapidly 
curved  around  into  the  spray  is  termed  the  spray  root.)  The  resulting 
surface  shape  expands  with  the  penetration  of  the  wedge  as  indicated  by 
the  several  positions  o'jajb.Cj,  o'2a2b?c^,  and  o'3a3b3c3. 

Under  the  assumed  perfect  properties  of  the  fluid,  such  a  series  of 
scalar  curves  for  a  constant  velocity  immersion  could  just  as  well  be 
represented  by  any  one  of  the  curves,  the  scale  of  the  diagram  then  de¬ 
pending  on  the  wedge  penetration.  Since  the  fluid  particles  are  not  re¬ 
stricted  in  size,  the  scale  of  the  diagram  is  of  no  consequence.  Thus, 
the  radiating  lines  oa,  ob,  etc,  define  corresponding  points  on  the  sur¬ 
faces  for  similitude  since  the  point  o  is  the  origin  of  the  motion. 

(It  is  interesting  to  note  that  several  points  other  than  the  origin 
chosen  could  be  used  as  a  center  of  symmetry  for  the  expanding  surface. 
The  wedge  apex,  y0,  the  spray  tip,  o',  or  the  intersection  of  the  wedge 
side  with  the  stilj  water  surface  could  be  considered  as  points  of  origin 
for  the  expansion  of  the  flow.  However,  the  latter  two  would  appear  as 
double  centers  for  the  complete  wedge,  and  all  three  havp  the  disadvantage 
that  the  center*  of  the  expansion  would  be  moving  with  respect  to  the 
initial  penetration  point  and  with  respect  to  the  fluid  at  a  great  dis¬ 
tance  from  the  wedge.  Since  such  motions  ox  the  center  of  expansion  would 
complicate  the  problem  unnecessarily,  the  initial  contact  point,  o,  will 
be  used . ) 
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Corresponding  points  on  the  surface  are  determined  by  the  ratios: 


0  laI  _  °' 2a2 

oaj  oa2 


Q't**  - 
ob2 


£3^ 

ob3 


e  tc. 


(13) 


Also,  the  velocity  at  corresponding  points  must  be  equal: 


kbV)  etc.. 


(14) 


where  ka,  kb,  etc  represent  the  proportionality  factors  between  the  velo¬ 
city  U  at  positions  a,  b,  etc.,  and  the  penetration  velocity, 

At  any  point  on  the  surface,  say  c3,  the  velocity  a*  determined 
from  consideration  of  this  condition  of  similarity.  From  the  instant  of 
contact,  the  surface  may  be  considered  to  have  been  sliding  past  the 
radius  oc  in  a  direction  tangent  to  the  surface  at  c3,  at  a  rate  suffi¬ 
cient  to  maintain  the  ratios  in  Equation  (13).  Thus,  a  particle  at  the 
surface  (at  point  c 3)  is  momentarily  moving  with  the  expansion  of  the 
radius,  and  at  the  same  time  sliding  past  the  radius  towards  the  spray 
region.  The  vector  sum  of  these  two  motions  (t/r?  +  USe)  yields  the  re¬ 
sultant  velocity  Uc  of  the  particle  at  the  point  m  question. 

The  radial  expansion  velocity,  Ur ,  and  the  surface  velocity  past 
the  radius,  Us ,  as  well  as  the  resultant  velocity,  U,  are  each  directly 
proportional  to  the  penetration  velocity,  V  (Equation (14)).  Thus,  for  con¬ 
stant  penetration  velocity,  V,  the  velocity  components  {Ure  and  Us.)  at 
position  c  are  constant,  and  the  vector  sum  of  the  velocities  may  be  ob¬ 
tained  from  the  radial  distance  and  the  arc  length: 


The  radial  distance  oc3  =  r  is  giver,  by: 


The  vector  s  is  constructed  tangent  to  the  surface  at  c  so  that  s 
equals  the  arc  length  o'c: 


The  vector  sum  is  then. 

s  +  r  -J*  Vc  dt  *  Uct  ,  (15) 

o 


Now  the  actual  fluio  particles  on  the  surface  do  not  follow  the 
radiating  lines  but  remain  at  a  fixed  arc  length  from  the  split  edge  of 
the  surface,  o'.  For  example,  a  particle  initially  at  rest  at 
rises  first  as  a  part  of  the  wave  surface  and  later  becomes  a  part  of  the 
spray.  The  velocity  of  the  particle  varies  with  its  position  on  the  sur- 
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face.  The  motion  of  the  particle  is  thus  a  function  of  its  relative  posi¬ 
tion  and  is  given  at  any  time  by  the  integral  relation: 

t 

V)  =  (€i  v)  dt  ,  (16) 

or  more  specifically, 
t 

*  "  ^  ~~J'u(f,v)  dl  (16a) 

l 

V  ~fv(Z.n)  dt  *  (16b) 

The  problem  is  to  determine  the  velocity  of  the  particle  (£,17)  at 
any  point  in  its  path,  and  then  to  relate  the  position  and  time  so  that 
he  integral  may  be  solved.  This  does  not  appear  to  have  an  analytic  so¬ 
lution,  but  the  graphical  solution  is  easily  carried  out. 

The  velocity  of  the  particle  {$' ,  tj)  may  be  found  directly  at  any 
particular  time  by  constructing  the  surface  of  the  proper  penetration 
locating  the  point  according  to  its  arc  length  from  o',  and  performing 
he  vector  addition  of  Equation  (15).  However,  this  would  be  a  tedious 
process  to  determine  enough  points  for  step-by-step  integration  of  Equa¬ 
tion  (16).  A  more  simple  alternate  method  is  to  change  the  variable  from 
time  to  arc  length,  which  is  done  in  the  following  manner. 

Assume  particle  (£,7?)  at  time  t2  to  be  at  position  c?.  Then,  at  t2, 
Ur:  '  Uc.v  whlch  may  be  determined  on  the  ts  curve  at  c3  .  In  other  wordq 
the  surface  points  djc3b3  correspond  to  earlier  times’  on  the  particle 
path.  The  time  t2  is  related  to  t3  by  the  arc  lengths  at  corresponding 
points  c 2  and  c3’  b 


t 


2 


ur,  in  general,  with  t  =  the  earlier  time  t,,  and  s 
arc  length  on  the  surface  at  time  t 3: 


the  corresponding 


and 


j 


& 


2 


dt 


ds 
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Substitution  for  dl  in  Equations  (16a)  and  (16b)  gives 


(17a) 


(17b) 


The  horizontal  and  vertical  components,  ut  and  vt ,  may  be  found  at 
any  point  along  the  surface  by  using  Equation  (15),  and  a  step-by-step 
integration  of  the  above  formulae  may  be  accomplished.  Any  trial  curve 
which  satisfies  the  continuity  condition  may  be  used  as  a  starting  point 
for  the  determination  of  the  surface  velocities.  The  solution  of  Equations 
(17a)  and  (17b)  determines  a  new  curve  which  is  used  in  turn  bo  obtain 
revised  valued  of  the  velocity  components.  This  iterative  process  is  coni 
tinued  until  negligible  differences  occur  between  computed -and  trial 
curves  of  surface  shape. 


Actually,  the  integral  equations  relate  the  rise  at  a  particular 
point  to  all  the  points  more  distant  from  the  wedge.  Thus,  the  path  of 
each  point  on  the  surface  follows  the  same  related  pattern,  and  the  sur¬ 
face  shape  fits  the  requirements  for  similarity.  Although  rigorous  ana¬ 
lytical  proof  is  omitted,  it  appears  from  the  actual  computations  which 
have  been  made  that  the  integral  relationship  (Equations  (17a)  and  (17b)) 
ggduces  the  number  of  possible  spray-wave  combinations  to  one  for  each 
spray  thickness. 


POTENTIAL  and  STREW  FUNCTION  ON  THE  FREE  SURFACE 


For  use  in  the  application  of  the  third  condition,  which  follows, 
a  potential  function  that  satisfies  the  first  two  conditions  is  assumed 
to  exist  for  the  fluid  field  of  each  spray-wave  combination  (even  though 
such  an  assumption  can  be  true  only  for  the  one  free  surface  shape  which 
meets  all  the  conditions  of  the  problem).  The  values  of  the  potential 
and  the  stream  function  'P  on  the  free  surface  may  be  found  as  follows. 

In  Figure  3,  the  tangential  velocity,  ’^p/'bs,  along  the  surface  at  a 
position  c  is  given  by  -0L  +  Vre  cos£.  In  terms  of  the  potential  •/>,  the 

tangential  velocity  =  For  points  alone  the  surface  (at  fixed  im¬ 

mersion),  cos£  is  given  by  dr/ds.  (The  values  of  dr/ds  for  the  particular 
surface,  such  as  3,  should  not  be  confused  with  dr/ds  -  r/s  for  a  constant 
position,  such  as  c.)  Then 


5s 


cos  C  -  V 


dr 

ds 


(18) 


The  general  differential  equation  is  obtained  by  multiplying  Equation  (18) 
by  t  ds  and  dropping  the  subscript  c  (noting  that  tl's  -  s  and  tVr  r) 
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- 1  d  ds  ~  s  ds  -  r  dr 
Integration  gives 

-  vt  fsLzjd)  +  .4 
2 

Since  -  0  at  s  =  r  -  ®,  A  -  0  and 


(19) 

(20) 


<P  r  r 2  -  s 7  _  My  ~Rs 


The  stre-'.m  function  'p  on  the  free  surface  is  found  in  a  similar 
manner.  At  position  c,  the  normal  velocity  is  given  by  Ure  sin?.  In  terms 

of  ■■■!, 

-  -I t*  Ur  sin?  . 

OS  c 

Multiplication  by  t  ds  gives: 
t  rls  -  t  sin?  ds  . 

The  integration  of  the  right-hand  side  is  accomplished  graphically  since 
it  represents  the  area  swept  out  by  the  radius  r.  Thus, 


,  rc 

J  t  jik.  ds  ~  -2  (prea  oo'co) 


and  in  general  terms,  the  general  position  s  replaces  the  particular  point 
c ,  and 


t  'Ps  t  'PQt  -  2  (area  oo'so) 


For  vi'  0  on  the  centerline  of  symmetry,  tip0  is  determined  by  the 
constant  normal  component  of  the  velocity  at  the  wedge  surface,  Un : 

Vn  1  V  cos/?  , 

t  %’  -  1  (>'0°')  ^  cos/5  •=•  (yQo')  y0  cos P  s  2  (area  oo'yQo) 
and  thus 


f  'Ps  -  2  (area  oo'yQo  -  area  oo'so)  .  (22) 

It  will  be  noted  that  as  s  -  «,  the  area  oo'so  becomes  equal  to  the 
area  oo'y0o  by  virtue  of  their  common  area  and  the  condition  of  continu- 
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ity.  Thus,  <As  =  0  =  4>w  ,  which  is  as  it  should  be. 

J  o 

I RROTATI QUALITY 


Irrotational  potential  flow  is  the  final  condition  for  the  field  of 
the  wedge  penetrating  the  fluid  surface.  Everywhere  within  the  fluid  and 
on  the  boundaries,  the  requirement  must  be  satisfied  that  V2<£  =  0.  This 
may  be  written  in  several  equivalent  forms: 


OX 


b2<p  b2ef)  +  b2ep  c‘Un 

by2  dn2  bs 2  bn  ds 


(23) 


From  the  graphical  constructions  made  for  the  condition  of  simi¬ 
larity,  it  is  possible  to  obtain  4>,  Un,  Us,  and  even  b(Js/bs  along  the 
free  surface  boundary;  bu  t  bUj  dr.  is  not  obtainable,  and  no  interior 
values  may  be  found.  Therefore,  it  is  not  possible  to  use  Equation  (23) 
as  a  direct  check  on  the  surface  shapes  obtained  through  application  of 
the  first  two  conditions.  Instead,  an  indirect  method  of  approach  is  used, 
based  on  Greene’s  theorem  (Reference  9). 

Greene’s  theorem,  which  relates  surface  integrals  to  volume  inte¬ 
grals  for  superimposed  potential  flows,  states  that  for  irrotationalpd- 
tential  flows: 


(24) 


where  and  are  any  two  irrotational,  potential  flows  superimposed  in 
a  common  region  S,  and  the  integral  is  taken  over  the  boundary  surface  of 
the  region  only.  If  the  known  irrotational  potential  function  for  a 
doublet,  is  chosen  for  4>t,  and  <f>7  =  (the  questionable  value  deter¬ 
mined  from  the  free  surface  shape),  the  integrals  will  check  only  if  4>  is 
also  a  potential  function  for  an  irrotational  flow. 

The  doublet  may  be  located  anywhere  on  the  assumed  free  surface. 
However,  the  field  of  a  doublet  decreases  rapidly  with  distance,  and 
the  choice  of  location  for  this  checking  purpose  may  be  quite  critical. 
Since  the  largest  deviation  from  irrotational  potential  flow  could  be  ex¬ 
pected  in  the  spray-root  region,  the  doublet  is  located  in  this  region  as 
shown  in  Figure  4.  (It  is  necessary  to  add  the  image  doublet  above  the 
wedge  to  give  bf^fbn  =  0  or.  the  wedge  surface,  because  the  graphical  con¬ 
struction  which  satisfies  the  similarity  condition  does  not  yield  any  in¬ 
formation  on  the  potential  along  the  wedge.) 

Consider  the  flow  i  the  immediate  area  of  the  doublet  as  shown  in 
Figure  5,  an  enlarged  section  of  Figure  4.  The  field  of  the  doublet  is 
given  by 


sina 

H 


(25) 


(26) 


i mr 
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from  which  the  velocities  due  to  this  field  are  piven  by 

.  .  sina 
dlt  liz 

d  „  cosa 
ZlBa  ft  2 

At  the  free  surface,  the  normal  component  is  given  by 


(27) 


J  sin6-~  -  cos8~~  =  — (-sin6  sina  -  cos6  cosa) 
cn  dR  Rda 


cos  (a  -  0) 
K2 


(28) 


The  potential  and  velocity  equations  for  the  image  doublet  are  the  same 
as  the  above  except  that  a  end  6  are  interpreted  as  positive  in  the 
counter-clockwise  direction  (and  of  course  the  R  values  are  measured  to 
a  different  center). 


On  the  4,  below  the  wedge  apex,  the  doublets  do  not  exactly  balance 
because  of  the  change  m  angle  between  the  wedge  bottom  and  the  4-  Since, 
at  this  distance,  the  velocities  due  to  the  doublet  or  its  image  are  small 
(as  1/R2),  the  net  difference  normal  to  the  %  is  very  small  and  may  be 
safely  neglected. 


The  surface  integral  required  to  check  the  V2<i>  condition  may  be 
evaluated  from  the  known  potential  field  of  the  doublets  and  the  measured 
field  characteristics  of  the  assumed  surface  shape  by  breaking  the  sur¬ 
face  into  a  number  of  segments: 


XkO 

/=y  s-o  s-s.-t  s  ,+e  s=®  ‘fc® 


x-o 

y*y * 


The  evaluation  of  the  integral  on  each  segment  will  be  discussed  separately 
arid  reduced  to  its  simplest  form  before  writing  Equation  (24)  in  its  en- 
t i rety . 

The  portion  of  the  surface  integral  on  the  &o£  symmetry  below  the 
wedge  apex  is  represented  by 
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This  portion  of  the  integral  may  be  shown  to  be  convergent  to  a  small 
value,  which  may  be  safely  neglected.  Over  the  part  of  the  surface 
B$/Bn  =  0,  is  small  (on  the  order  of  l/l 0,  and 


J *  (0)  (sin a/R)  ds  c  0 


for  s  -  R  at  infinity.  Also,  for  one  doublet  is  very  small  at  and 
below  y0  (on  the  order  of  1/fi2).  For  the  pair  of  doublets,  the  value  is 
the  difference  between  the  two,  an  evei  smaller  value.  On  the  below 
the  wedge,  4>  is  not  given  by  the  analysis  so  far,  but  from  the  expanding- 
plate  analogy  (Reference  3),  the  value  is  on  the  order  of  Ry0  at  the  wedge 
apex,  and  diminishes  rapidly  below  the  wedge,  giving 


f  Rv  .£.?±S*-zJLds  , 
*  y°  R2 


which  is  convergent  for  s  —  R  -  °°.  Since  the  small  finite  value  which 
might  accrue  from  this  integral  appears  to  be  much  lower  than  the  mechan¬ 
ical  accuracy  of  the  graphical  method  to  be  used,  this  term  may  be  neg¬ 
lected. 

The  value  of  the  integral  along  the  solid  boundary  from  wedge  apex 
to  spray  tip. 


*  »  a 

may  be  determined  analytically.  Along  the  wedge,  <t>  is  not  known,  but  is 
finite,  while  B0^/B n  -  0,  so  that 


y  4> - ^  Js  =  0. 
y-  B  n 

*  O 

The  doublet  and  its  image  add  along  the  wedge: 

rt .  ~  n  sina 
*d  ■ 

Since  Bf>/'Bn  =  V  cos,£  and  ds  =  R  do./s ina, 
o'  i w' 

f* ■£  .  ,/ c  =  /*  9  V  rna/i  da  =  9  V  /,  roc£  . 

J  a  Bn  V  -  t  -  ■  , 

ay„ 

where  <u  is  yQdo', 

Two  terms 


sEs 

►  d 


which  represent  the  integrals  on  the  free  surface  from  the  spray  tip  to 


the  doublet  and  from  the  doublet  to  infinity  are  indicated  because  the 
free  surface  is  theoretically  broken  for  a  distance  2c  at  the  doublet 
Is  «  s^),  At  all  other  points  of  the  free  surface,  4>  and  are  given 
by  constructions  for  the  assumed  surface,  and  and  ctyj/'dn  are  given  by 
Equations  (25)  and  (28).  Because  of  the  discontinuity  at  the  doublet 
center,  a  small  semi-circular  path  is  followed  around  the  singular  point 
to  obtain 

V'- 


(The  distance  is  taken  so  small  that  the  effect  of  the  image  doublet 
is  infinitesimal.)  The  doublet  potential  and  normal  velocity  components 
are  given  by  Equations  (26)  and  (27).  The  potential  and  normal  velocity 
of  the  main  flow  field  are  obtained  on  this  small  semi-circular  segment 
by  assuming  no  var lation  of  the  field  velocity  for  the  small  area  involved. 
As  shown  m  Figure  5. 


'A  *  R  cos  (for  R  -  e  only) 
d  ™ 


and 


dfl 


COSO.  | 

tfn  ] 

d 


sino-|^  I 

os  I 
d 


where  |  ^  indicates  a  constant  value  for  the  main  flow  field  in  the  im¬ 
mediate  region  of  the  doublet.  Thus, 


4>  -■  <t>\  *  R  cos'*t2£l|  -  R  sina-Mj 

id  I  <>s  I 

d  d 


and 


M  ds 

11  r)n 


-tt/2 

•/ 


-/  (^f) 

d  d 

.  v 
"  2  Ss 

d 


Also; 
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a 

f  ■ 


~bt, 

'P-~  ds 


"l 


'/( 


4>  +  B  cosa-M|  -  R  sina-|^| 

d  on  I  as  I 

d  d 


_3#|  \  /  -sinaV 
3s  1  A  R2  I 


-  +  JL 

2  3s 


The  closing  segment  of  the  contour  from  the  free  surface  to  the  <£, 
of  symmetry  on  the  infinite  boundary  is  represented  by 

/ 

Even  though  this  integral  is  over  an  infinite  boundary,  it  is  definitely 
zero  because  <P,  fy,  3<p/3n,  and  3^>d/3n  are  zero  on  the  infinite  boundary. 

Now  that  the  integral  for  'each  segment  of  the  contour  has  been  re¬ 
viewed,  it  is  simple  to  write  the  integral  Equation  (24) 


A-t?* 


which  may  be  broken  down  into  sections  according  to  the  discussion  above 
to  give; 

’/  •  <3*» 


This  reduces  to- 


*  2  V  n  cos/5  - 


CD 

A 


ds  . 


Equation  (31)  expresses  the  condition  of  irrotationality  in  a  form 
which  may  be  applied  to  the  trial  curves  which  have  been  obtained  from 
the  first  two  conditions.  By  trial-and-error  solution,  the  final  shape  of 
the  fluid  surface  which  will  satisfy  the  third  and  final  condition  may  be 
found  for  each  wedge  angle. 

DETERMINATION  OF  LOADS  and  FftESSURES 

From  the  above  portion  of  the  analysis,  the  shape  of  the  free  sur¬ 
face  for  potential  flow  due  to  the  constant-velocity  immersion  of  a  wedge 
has  been  determined.  The  values  of  velocity,  potential,  and  stream  func¬ 
tion  on  the  free  surface  corresponding  to  this  shape  are  directly  ob- 
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tamed,  as  indicated  previously.  However,  (in  order  to  obtain  the  pres¬ 
sure  on  the  wedge)  velocity  and  potential  values  along  the  wedge  surface 
are  required,  and  these  are  not  directly  given  by  the  analysis  so  far. 
The  study  will  be  extended  to  yield  this  information,  but  first,  the 
method  of  obtaining  the  total  load  as  proposed  by  Wagner  (Reference  8) 
will  be  developed. 

WAGNER'S  METHOD  FOR  TOTAL  LOAD 

Penetration  of  the  wedge  into  the  fluid  generates  pressures  on  the 
wedge  which  appear  in  the  fluid  flow  as  momentum  and  kinetic  energy  (under 
the  initial  assumptions  of  the  perfect,  fluid  without  restraining  bounda¬ 
ries).  The  total  force  at  any  time  equals  the  rate  of  change  of  momentum, 
and  the  work  done  by  the  force  is  equivalent  to  the  kinetic  energy  of  the 
f 1 uid . 

The  momentum  of  the  fluid,  M,  is  given  by  summation  of  the  products 
of  incremental  volume,  d(vol),  and  the  local  velocity  vector,  V<£  =  grad<£: 

■V  ^  £^v°D 

Since  the  local  velocity,  V<£,  at  any  position  is  constant  for  a  constant 
velocity  penetration,  the  change  in  momentum  with  respect  to  time  at  a 
given  position  in  the  fluid  is  given  by  the  rate  of  increase  of  the  dif¬ 
ferential  volume.  For  the  two-dimensional  case  under  consideration,  the 
differential  volume,  d(vol),  to  occupy  a  given  space  increases  directly 
as  the  area  of  immersion,  or  directly  as  the  square  of  the  depth  of  im¬ 
mersion.  Thus,  d(vol)  Z  A  Z  yQ2  z  t2  . 

The  total  force  is  then  given  by 


where  V<£  and  J(vol)/f2  are  for  relative  positions  and  independent  of  time 
t.  The  differential  of  the  integral  with  respect  to  time  gives 


P 


d(vol)  . 


(32) 


This  integral  can  be  transformed  to  a  surface  integral  by  the  theorem  of 
divergence  (Gauss’  theorem) 


J(vol) 


where  the  differential  area,  df,  is  a  vector  normal  to  the  surface,  posi¬ 
tive  inward  with  respect  to  the  surface.  Substitution  in  Equation  (32) 
g  i  ves 
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Pt  *  4>  df  ,  (33) 

which  for  the  two-dimensional  case  under  consideration  may  be  written: 

Pt  =  -2pJ*<}>  ds  , 

where  P  represents  the  load  per  unit  length  of  the  infinite  wedge. 

Because  of  the  symmetry  of  the  problem,  the  horizontal  components 
of  the  vector,  df,  balance.  Thus,  only  the  vertical  components  need  be 
integrated  to  obtain  the  load  P  (which  must  obviously  be  normal  to  the 
initial  fluid  surface  for  the  symmetrical  problem  under  consideration). 
The  vertical  components  of  the  incremental  vectors,  <b  ds ,  may  be  written 
as  4>  dx,  where  dx  is  the  horizontal  projection  of  an  elemental  surface 
length.  The  direction  (up  or  down)  of  4>  dx  is  determined  by  the  sign  of 
the  quantity  dx,  depending  upon  the  position  of  the  surface  with  respect 
to  the  fluid,  as  shown  in  Figure  6.  Along  the  wedge,  dx  =  cosp  ds  and  the 
vertical  component  of  4>  ds  =  cos/3  ds.  On  the  spray  surface  the  inward 
vector  is  upwardly  inclined  so  that  ds  ~  -dx ,  while  on  the  wave  surface 
the  inward  normal  is  downward  so  that  ds  %  dx . 


Now,  if  Vn  represents  the  component  of  the  field  velocity,  V ,  normal 
to  the  surface  and  positive  inwards,  Vn  will  be  negative  in  the  spray  and 
positive  in  the  wave  area  as  shown  in  Figure  6.  Thus,  in  both  regions,  dx 
may  be  written  as 


and  the  vertical  component  of  df  is 
4>  dx  -  ds 

The  integral  for  total  load,  Equation  (33),  may  then  be  written  in 
two  parts:  one  for  the  body,  Sg,  and  one  for  the  free  surface,  5^; 


-2/?y~  $  ds  *  ^ 


A  >•  K. 

Pt  =  2p  cospj  <P  ds  -'IpJ  d>y  ds  .  (341 

.  ** 

The  kinetic  energy,  K.E.,  of  the  fluid  field  is  given  by  Greene’s 
theorem  in  terms  -f  a  surface  integral  (Reference  9): 


■  °/sM 


+ - 2-  )  rf(vol)  3  --£■ 


where  positive  energy  is  indicated  by  flow  outward  through  the  surface, 
The  kinetic  energy  per  unit  thickness  cf  the  wedge  is  therefore 
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K.E.  -  -4  ds  . 

2%/  dn 


By  virtue  of  the  uniform  growth  of  the  loaded  area  for  the  simple 
wedge  immersing  at  constant  velocity,  the  load  on  the  wedge  increases 
linearly  with  time.  (This  may  also  be  deduced,  mathematically,  from  Equa¬ 
tion  (32).)  Thus  P  -  (P/t)t,  and  the  work  done  is  given  by  the  integral 
of  the  load,  P,  times  the  distance,  dy  (  =  V  dt): 

Work  ■-  f  P  dy  t  V  dt  s{~j  V  -y  . 

Since  the  work  done  on  the  perfect  fluid  must  equal  the  change  in  kinetic 
energy, 

ill  -  K.E.  -  foMds  .  (36) 

2  2  J  on 

s 

Since,  at  the  body,  ^4>/2n  =  +  V  cos/3,  and  on  the  free  surface,  o<pj'0n  is 
negative,  variable  with  position,  Equation  (36)  may  be  written: 

~f  (+  V  c°s ef  *ds  fds) 

% 

or 


Pt 


p  cosjS J'  4>  ds 


ds 


(37) 


It  will  be  noted  that  both  Equations  (34)  and  (37)  contain  a  term 
which  requires  values  of  the  potential  on  the  body  for  solution.  However, 
they  may  be  solved  simultaneously  to  yield  a  solution  for  the  Load  P  in 
terms  of  the  free  surface  integrals  alone,  eliminating  the  troublesome 
term. 


Equation  (34)  minus  twice  Equation  (37)  gives: 
Pt  -  2  Pt  -  -2 p J*  <t>~  ds 


- PVt  *  2  c 


As[*j*tds  's[,pv" ds) 


-PVt  =  2 0 


(38) 


This  integral  may  be  evaluated  for  each  deaarise  investigated  since 
the  required  potential  and  normal  velocities  are  given  at  every  point  on 
the  free  surface.  In  general,  the  velocity  and  potential  have  been  ob- 
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tained  in  terms  of  distances  in  the  geometric  construction  ( tV ,  t<p) . 
Equation  (38)  may  be  converted  to  this  form  by  multiplication  by  t2: 


-2p /  t<A 


n  dn 


or  since  y0  = 


Equation  (39)  is  the  final  form  for  convenient  evaluation  of  the 
load  on  the  immersing  wedge  in  terms  of  the  free  surface  shape  determined 
in  the  previous  section. 

It  may  be  seen  by  reference  to  Equation  (20)  and  Figure  3  that  the 
major  portion  of  this  integral  accrues  in  the  spray  region.  The  value  of 
potential  is  large  at  the  spray  tip  and  decreases  rapidly  to  low  values 
past  the  spray  root,  while  the  relative  normal  velocity  component  in¬ 
creases  to  a  much  smaller  extent  in  the  spray-root  region.  The  product  of 
potential  and  normal  velocity  diminishes  rapidly  to  zero  along  the  wave 
surface  towards  infinity  because,  while  the  velocity  approaches  a  con¬ 
stant  value  [iV'n  -  t(30/3n,)  =  iV],  the  product  t4>  diminishes  as  the  dis¬ 
tance  squared . 

Closer  scrutiny  of  the  relative  velocity  [Vn  -  (d<p/on))  shows  the 
direct  relationship  between  the  spray  thickness  and  the  total  force  ac¬ 
cording  to  Equation  (38).  The  velocity  [ Vn  -  (30/3 n)3  is  actually  the  nor¬ 
mal  component  of  the  vector  sum  of  V  and  Uc  in  Figure  3.  In  the  spray 
region,  this  is  practically  constant  and  equal  to  the  velocity  of  the 
spray  tip  along  the  wedge  times  the  angle  between  spray  and  wedge.  Thus, 
the  thickness  of  spray  determines  the  relative  normal  component  in  the 
region  in  which  the  greatest  part  of  the  force  integral  is  determined, 
and  so  the  force  is  roughly  proportional  to  the  spray  thickness. 

DERIVATION  OF  POTENTIAL  AND  VELOCITY  DISTRIBUTIONS  ON  THE  WEDGE  SURFACE 

From  the  simple  graphical  constructions  described  in  the  previous 
section,  the  analysis  of  the  free  surface  shape  for  any  wedge  angle  will 
give  directly  the  velocity  and  potential  distributions  on  the  free  sur¬ 
face.  Also,  a  method  has  been  given  for  obtaining  the  total  load  on  the 
wedge,  which  circumvents  thenecessity  for  detailed  flow  information  along 
the  wedge  itself.  However,  the  distribution  of  the  pressure  on  the  wedge 
has  not  been  given,  and  for  that  purpose  the  potential  and  velocity  dis¬ 
tributions  on  the  wedge  surface  are  required. 

On  the  premise  that  the  uniqi  -»  free  surface  rhape  for  each  wedge  is 
a  result  of  potential  flow,  the  velocity  and  potential  at  any  point  may 
be  represented  by  complex  analytic  functions: 


to  which  the  Cauchy  integral  equation  may  be  applied,  relating  the  point 
value  of  a  function  to  the  integral  around  a  bounding  contour  of  values 
of  the  function  at  other  points.  The  standard  integral  for  the  point 
within  the  boundary  is' 


f(*o>  = 


±-[jLL 

■vj  2  ’ 


ii—  dz  , 


where  z  is  the  complex  coordinate  of  the  position  |n  the  field.  Since 
points  on  the  solid  boundary  are  to  be  investigated,  the  bounding  contour 
circles  only  halfway  around  the  point  in  passing  so  that  the  factor  is 
reduced  to  v,  and  the  principal  value  of  the  integral  is  used.  Thus,  for 
the  velocity  U  =  tip  -iv p  *>  f(z): 


0<zo) 


ViJl  ” 


lL  dz 


This  equation  may  be  separated  into  real  and  imaginary  parts  by 
placing  the  origin  of  coordinates,  Zp  -  xp  +  iy  p,  at  the  point  zQ  being 
investigated,  and  writing  z  in  the  polar  form: 


2  =  pei6  *  eln P*i0 


dz  =  ew^ie  [d  fin  p)  +  id  6]  =  2  [d(ln  p)  *  i  d  0]  . 


Substitution  of  these  expressions  for  z  and  dz  into  Equation  (41)  gives 


U(*o>  ~  uop  -  lv 


Rearranging  the  terms  gives 


°R  vi  J 


[tip  -  iv q]  [dfln  p)  +  i  d  6] 


iv  ar.  +  mn 
°0  °/3 


f 


•dp  d(lnp)  *  Vp  d&  +  1  [u&d6  •  Vp  dfln  4)3  , 


;o  that  the  real  and  imaginary  parts  of  U(r-a)  may  be  separated: 


%  ‘■i/'* df>  *  d( In  p) 
...  -ir*.  dO  -  Vn  cffln  p) 


(45a) 


(45b) 
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The  construction  for  performing  this  integration  is  shown  in  Figure 
7.  The  complex  coordinates  are  taken  with  the  iyc  axis  along  the  wedge, 
positive  towards  the  spray  tip.  The  Xp  axis  is  then  normal  to  the  wedge 
surface,  positive  into  the  fluid.  The  axes  are  so  chosen  that  the  integral 
for  V0  need  not  be  evaluated  along  one  side  of  the  wedge  since  u.p  is  con¬ 
stant  (=  V  cosfi)  and  d6  =  0  (except  for  the  infinitesimal  region  of  the 
singularity).  On  the  far  side  of  the  wedge,  the  value  of  Vp  must  be  esti¬ 
mated  until  the  first  solutions  for  vp  are  obtained.  However,  these  points 
are  sufficiently  distant  to  have  small  effect,  and  a  rapidly  convergent 
solution  is  obtained. 

The  potential  on  the  wedge  surface  may  be  obtained  by  integrating 
the  velocity  along  the  wedge  from  the  spray  tip  (where  the  free  surface 
formula,  Equation  (20),  gives  the  coincident  value)  to  the  apex  of  the 
wedge.  Tne  potential  may  also  be  found  directly  by  the  integral  relation 
corresponding  to  Equation  (45)  (<f>  in  place  of  up,  and  </»  in  place  of  -Vp). 
This  is  more  laborious,  and  only  a  few  points  need  be  investigated  to 
check  tiie  results  from  the  velocity  integration  along  the  surface  of  the 
wedge. 

PRESSURE  DISTRIBUTION 

The  velocity  distribution  and  potential  along  both  the  free  surface 
and  the  immersed  portion  of  the  wedge  have  been  determined  in  the  previ¬ 
ous  sections.  It  now  remains  to  compute  the  pressures  which  correspond  to 
this  potential  flow. 

As  was  noted  before,  particles  in  the  fluid  move  relative  to  the 
shape  of  the  flow  field  (across  instantaneous  streamlines).  Thus,  the 
simple  Bernoulli  pressure  equation  may  not  be  u  .ed  since  particular 
streamlines  exist  only  instantaneously,  even  though  the  flow  pattern  is 
always  similar.  Fcr  this  problem,  the  general  equation  for  fluid  pressure 
in  a  potential  flow  muot  be  used  (Reference  9): 

JL  -M-Ol  <  F(t)  +  Cl  .  (46) 

P  2 

Since  che  fluid  has  been  assumed  to  be  at  rest  at  infinity,  with  gravity 
neglected,  the  gravitational  field  and  extraneous  impulses,  represented 
by  Cl  and  F(t),  ere  zero.  The  resultant  field  velocity,  V,  has  been  deter¬ 
mined  on  the  free  surface,  (Equation  (15)),  and  on  the  wedge  (Equation 
(45)).  The  rate  of  change  of  potential  with  time,  ap/’dt ,  maybe  determined 
as  foil  n»'3  0 

At  any  position  r/yQ  in  the  fluid  (see  Figure  8),  the  potential  <£ 
is  a  linear  function  of  time  (since  time  is  directly  related  to  the  scale, 
or  size,  of  the  symmetrical,  constant-velocity  penetration  under  consid¬ 
eration)  r 


£ 


(47) 
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However,  at  a  fixed  point  (x,y)  in  the  field,  this  is  not  so.  Con¬ 
sider  the  point  (x,y)  to  coincide  w:'th  the  position  r^/y0  * j  •  At 
t  =  t{  +  At ,  the  position  rt/y0  is  a  ’istance  Ar  =  rAt/t  from  (x,y).  (For 
small  At,  the  velocity  U  in  the  region  may  be  assumed  constant  in  magni¬ 
tude  and  direction.) 

Then,  at  t  “  tJf 


at  t  «  ti  +  A t,. 


where  !J  cosa  is  the  component  of  the  velocity  along  r .  Thus, 

**  "  ***■*>, ,*tt '  *<*■>>,,  “  "  *r  costt "  ■ 

With  i  =  tj  *  At,  frorr.  liquation  (47), 

‘  *rt  ’*r,  • 

t  tj  t,  tj  l} 

Equation  (48)  then  reduces  to 


fyfr  r  +  (jy  ■  -  U  Ar  cosa 

t.  rt. 


t 


The  partial  derivative  ~&p}'dt.  is  obtained  from  the  limit  of  Equation 
(49)  as  At  approaches  0: 


~~  «■  lim  -M  lim  (  + - —  *-x  ”  cosa>) 

dt  At-o  At  At-o  \  tj  At  / 


w  *  jt.  -_L  y  cos#  f  (50; 

t  £ 

where  the  subscripts  are  no  longer  needed  since  the  analysis  applies  to 
all  points  in  the  fluid.  The  pressure  distribution  is  then  determined  by 


_H  -  -  jL  *  JL  V  cosa 
i o  t  t 
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where  the  expression  for  cx£/3t  from  Equation  (50)  has  been  substituted  in 
Equation  (46)  (with  fl  =  F(t)  =  0). 

On  the  free  surface,  the  pressure  should  be  zero.  This  may  be 
checked  by  referring  back  to  the  geometric  relations  between  the  velo¬ 
cities  and  potential  shown  in  Figure  3  and  Equation  (21).  These  values 
substituted  in  Equation  (53)  yield 


+  U  Ur 


cosa 


Ml 

2 


t 


the  right-hand  side  of  which  will  be  recognized  to  be  zero  from  the  solu¬ 
tion  for  a  triangle  side  Us  in  terms  of  Ur,  U,  and  the  included  angle,  a. 

On  the  wedge,  the  pressure  may  be  found  by  substituting  the  values 
of  r,  4>,  V,  and  a  determined  previously.  The  total  load  on  the  wedge  is 
obtained  from  the  integration  of  the  -edge  pressure  across  the  projected 
wetted  width.  This  value  may  then  be  compared  with  that  obtained  on  the 
basis  of  Wagner’s  analysis  for  the  integral  over  the  free  surface  only 
(Equation  (39)). 

Maximum-pressure  location  and  magnitude  along  the  wedge  may  be  found 
by  setting  equal  to  zero  the  first  derivative  with  respect  to  s  of  Equa¬ 
tion  (51): 


dP  _  1  dtp 

~Js  TUT 


*  —  U  cosa  Ml 
t  ds 


cos  a 


dU 

7s~ 


uM!L° 

ds 


0. 


Now  U  cos  a  =  d<j)/dr  so  that  at  P  m  (or  P  .  ), 

O  «  X  *  1  K 


—  cos  a.™  - 
t  ds  ds 


(52) 


A  maximum  or  minimum  can  occur  when  dU/ds  =  0  or  when  r/t  =  t//cosa.  Since 
the  former  occurs  at  the  spray  tip,  a  minimum  pressure  point,  it  may  be 
ignored.  Substitution  of  the  latter  into  Equation  (51)  gives 


P  «  _  _  jj.  +  11  ~ 
P  t  2 


(53) 


or  in  terms  of  r, 


P**x  c  __0  ^  ’•2  cos2*  (54) 

P  t  2t2 


Further  study  of  the  relationship  r/t  »  fj/ cosa  shows  that  at  this 
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point,  a  must  be  zero.  This  follows  from  the  observation  that  the  com¬ 
ponent  of  U  at  the  wedge  normal  to  the  wedge  is  constant  and  equal  to 
the  component  of  r/f  normal  to  the  wedge.  Thus,  r/t  can  equal  U/cosa. 
only  when  r/t  ~  U  and  a  °  0°.  Dividing  both  sides  of  Equation  (54)  by 
K2/2  =  y2/2t2  and  putting  cosa  =  1,  gives 


P 


■  AX 


P\'2/2 


at  the  point  (or  points)  where  r/t  =  U. 


(55) 


It  is  interesting  to  note  that  the  peak  pressure  point  r/t  =  U 
corresponds  to  a  stagnation  point  in  that  the  same  particle  occupies  the 
same  relative  position  at  all  times.  However,  the  pressure  is  not  the 
simple  relation  pU2/ 2,  but  is  greater  or  less  depending  on  whether  the 
potential  at  the  point  is  less  or  greater  than  0. 


RESULTS  OF  ANALYSIS  AND  COMPARISONS  WITH  OTHER  METHODS 

SURFACE  SHAPE 

The  details  of  computation  and  drafting  required  for  the  application 
of  the  continuity,  similarity,  and  irrotationality  conditions  to  the 
penetration  of  a  40°  deadrise  angle  wedge  are  given^in  Appendix  I.  Simi¬ 
lar  calculations  were  performed  for  angles  of  20°  ,  30°,  and  50°.  The 
final  surface  shape  determined  by  the  iterative  solution  is  shown  for 
each  angle  in  Figure  9. 

The  spray  thickness  appears  to  be  the  most  critical  characteristic 
for  a  quantitative  description  of  each  surface.  However,  since  the  spray 
is  tapered  to  a  point  at  its  edge  and  merges  into  the  spray-root  region, 
it  is  necessary  to  establish  a  consistent  definition  for  the  spray  thick¬ 
ness.  It  was  noted  that  in  every  case  the  spray  root  comes  quite  close  to 
the  nominal  half-width,  c.  Therefore,  this  nominal  width  was  used  as  a 
reference  line,  and  the  spray  thickness  was  defined  as  the  distance  from 
the  wedge  surface  to  the  intersection  of  the  straight  portion  of  the  spray 
surface  extended  to  the  half-width,  c,  as  shown  in  Figure  9.  The  ratio  of 
spray  thickness  to  half-width,  S/c,  defines  the  relative  spray  for  each 
wedge  angle.  The  three  surfaces  computed  in  detail  were  used  as  the  basis 
for  the  plot  of  relative  spray  thickness  vs.  wedge  angle  given  in  Figure 
10.  The  points  are  approximately  fitted  by  the  formula 

JL  --  .LiP/1  -tan£  .  (56) 

c  3 . 3n 


*  Tift  »pr»y  th,ckne«»  for  .9  -  20°  »»*  taken  fr  the  carve  determned  by  the  other  three,  end 
enly  eoatlauit,  end  iieilintf  condition*  *ere  eppiied. 


AOs 


In  References  3  and  8,  an  analytical  solution  is  obtained  for  the 
flow  in  the  region  of  an  immersing  wedge  by  comparison  with  the  solution 
for  the  flow  around  an  expanding  flat  plate.  This  solution  gives  an  ex¬ 
pression  for  the  rise  of  the  free  water  surface  (Equation  (3)),  with  no 
spray. 


The  mathematical  analysis  of  the  spray-root  region  (References  2 
and  8)  describes  the  flow  in  the  immediate  area  of  the  spray  root,  but 
the  boundaries  of  the  problem  are  doubly  infinite  and  do  not  correspond 
to  the  conditions  of  the  penetrating  wedge.  The  flow  is  defined  in  terms 
of  the  spray  thickness  which  is  not  readily  determined  for  the  finite 
problem. 

Wagner  obtained  an  approximate  expression  for  the  spray  thickness 
by  comparing  the  flow  fields  at  infinity  for  the  spray-root  and  expanding- 
plate  analyses  (Reference  8,  Equation  47): 


5  _ 


c 


where  u  =  2  tanfi/v,  and  therefore, 

1  *  tan-J  .  (57) 

c 


However,  the  logical  procedure  by  which  this  formula  was  obtained  is  not 
readily  apparent  in  the  referenced  report.  Later,  unpublished  work  by 
Wagner  indicated  that  a  factor  o£  0.7  (approximately)  should  be  applied 
to  the  spray  thickness  obtained  from  Equation  (57). 

A  more  reasonable  comparison  of  the  spray-root  and  exnanding-plate 
analyses  may  be  made  in  the  region  of  the  spray  root  itself.  This  com¬ 
parison  is  based  upon  the  assumption  thnt  the  stagnation  point  of  the 
spray-root  analysis  will  occur  at  the  peak -pressure  point  for  the  expand¬ 
ing-plate  analysis.  A  superposition  of  the  two  flows  with  the  'stagnation' 
as  a  common  point  is  sheen  in  Figure  11. 


In  the  previous  analysis  of  spray-root  flow  (Reference  2),  the 
scale  was  indeterminate  and  the  local  region  was  not  rented  t.o  the  over¬ 
all  picture  of  a  wedge-penetration  problem.  The  distance  A  in  Figure  11 
is  approximately  eq"al  to  the  distance  from  spray  root  to  stagnation 
point  of  the  previous  lysis,  and  may  be  obtained  in  terms  of  the  spray 
thickness,  — ,  from  ; igurc  5  of  Reference  2  by  projection  of  the  root  into 
the  *7t/§  axis  (at  a  value  of  6).  For  the  angles  investigated,  A  =  £w/S 
and 


c  -  x  -  ( S  tan#  +  A)  cos#  *  5  (ein.3  +  —  cos#)  .  (58) 

v 

The  expanding-plate  analogy  mentioned  previously  does  not  give  a 
practical  physical  solcticn  because  of  the  absence  of  spray,  but  by  that 
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method  of  analysis,  a  definite  relationship  is  established  between  the 
wedge  angle  and  the  stagnation  point.  Now  x/c  is  given  by  Equation  30  of 
Reference  3  as 

—  -  Jl  -  X2  ®  x/ i  -  (2  tunp/w)2  , 
c 

and  substitution  of  this  value  from  the  expanding-plate  analogy  into  Equa¬ 
tion  (58)  relates  the  spray  thickness  ratio  8/c  to  the  wedge  angle,  /3: 

s  -  c  sf  1  -  (2  tan/3/r.) 2  e  8  (sin/S  +  —  cosfi) 

77 


—  -  L  ~  ^  f  *•  (2  tanlAO2 
c  aivfi  +  i-  tosj3 


(59) 


Equations  (57)  and  (59)  are  plotted  in  Figure  10  for  comparison 
with  the  points  determined  by  the  present  analysis.  It  will  be  noted  that 
the  three  points  computed  in  detail  (at  30°,  40°,  and  50°  deadrise)  ap¬ 
proach  the  fcurve  of  Equation  (59)  at  the  lower  angles.  Also,  Equation 
(59)  approaches  70%  of  Equation  (57)  at  low  deadrise  angles. 


VELOCITY  POTENTIAL,  AkO  STREAM  FUNCTION 


Under  the  conditions  for  inotational  flow  of  a  perfect  fluid,  the 
shape  of  the  free  surface  hss  been  determined  for  several  wedge  angles. 
In  fact,  the  solution  for  the  surface  shape  is  the  key  to  the  field  of 
fluy  everywhere  in  the  fluid.  The  values  of  the  velocity,  potential, 
and  stream  function  on  the  free  surface  are  given  directly,  while  the 
values  on  the  wedge  must  be  estimated  and  then  refined  by  an  iterative 
solution. 

FRht  Sl'RFACf 

Velocity  and  potential  at  the  fret  surface  are  given  by  Equations 
(IS)  and  (21)  respectively.  These  have  already  been  computed  for  the 
several  stages  of  the  solution  for  the  free  surface-  shape.  The  final  set 
of  values  for  each  wedge  angle  are  plotted  in  Figures  12,  13,  and  14  as 
U/V,  a,  and  <£t/c7  respectively  vs,  position  on  the  surface,  s/c.  The 
stream  function  is  given  on  the  free  surface  by  Equation  (22).  The  values 
of  yn/c2  vs.  s/c  for  the  several  wedge  angles  are  plotted  in  Figure  15. 

SOL  1C  BOUNDARY 

At  the  wedge,  the  normal  component  of  the  velocity,  up,  is  constant 
(V  cos#),  and  the  stream  function  is  given  directly  by 
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<P  B  sV  cos/3  *  y^o'  V  cos/3, 

where  s,  the  distance  from  the  spray  tip,  will  be  negative  along  the 
wedge.  The  tangential  component  r.i  the  velocity  at  the  wedge,  vp,  must  be 
estimated  so  that  the  complete  integral,  Equation  (45),  may  be  evaluated. 
The  revised  values  obtained  from  the  iterated  solution  of  the  Cauchy  in¬ 
tegral  for  each  wedge  angle  are  plotted  in  Figure  16.  These  data  are  com¬ 
bined  with  the  constant  up  and  shown  in  Figures  12  and  13. 

Values  of  the  potential  on  the  wedge  were  obtained  from  the  inte¬ 
gration  of  the  velocity  distribution  from  wedge  apex  to  spray  tip  and 
vice  versa.  The  end  points  for  the  potential  (at  spray  tip  and  wedge  apex) 
were  obtained  from  Equation  (21)  and  the  potential  relationship  corre¬ 
sponding  to  Equation  (45).  c-ince  the  mechanical  work  involved  in  the  con¬ 
struction  practically  precludes  a  perfect  agreement,  the  curves  of  4>t/c2 
vs.  s/c  along  the  wedge  presented  in  Figure  14  represent  average  values 
obtained  by  working  from  either  end  toward  the  center  of  the  wedge  side. 

It  is  of  interest  to  compare  the  potential  and  velocity  from  the 
present  analysis  with  those  for  the  immersing  wedge  as  given  by  the  ex¬ 
panding-plate  analogy.  Equation  4  of  Re ference  3  gives  for  the  potential 
on  the  plate 

<p  =  - v/c 2  -  x2 


or 


*L-  -h.!i  -  (JL)2 

c  2  C  V  \  C  l 

This  expression  is  indicated  by  the  dotted  lines  in  Figure  14  (with  signs 
reversed  for  clarity). 

Equation  5  of  Reference  3  gives  for  tne  velocity  along  the  plate 

v,  --  U  =  V,  x 

- - 

v/c2  -  X2 
or 


J(c7x)2  -  1 

This  expression  is  represented  by  the  dotted  lines  in  Figure  16. 

It  will  be  noted  that  the  velocity  and  potential  distributions 
along  the  wedge  obtained  by  the  earlier  method  (Reference  3)  are  in  par¬ 
tial  agreement  at  the  lower  angles,  but  that  the  results  of  the  twometh- 
ods  diverge  at  higher  angles. 
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PRESSURES  AND  LOADS 


The  pressure  distribution  on  the  wedge  surface  is  given  by  Equation 
(51).  The  values  of  <j>t/c7,  U/V,  and  a  have  been  calculated  already  and 
are  plotted  in  Figures  14,  12,  and  13,  while  r/Vt  may  be  measured  direct¬ 
ly.  The  computed  pressures  for  the  20°,  30° ,  40°,  and  50°  deadrise  wedges 
are  plotted  in  Figure  17. 


Total  load  on  the  wedge  may  be  obtained  by  integrating  the  pressure 
distribution  along  the  bottom.  Because  of  the  symmetry  of  the  wedge,  only 
the  vertical  component  need  be  taken  since  the  lateral  forces  from  the 
two  sides  will  cancel  each  other.  Thus,  for  the  complete  wedge  (two 
sides) : 


and 


P 


2  c  os  A 


ds 


The  loads  obtained  by  integrating  the  pressures  for  the  three  wedges  in¬ 
vestigated  have  been  used  to  construct  the  curve  of  P/c(pV3/ 2)  vs.  A 
given  in  Figure  18. 


For  comparison  with  Wagner’s  analysis  of  the  total  load,  the  values 
of  P/c(pV2/ 2)  obtained  by  Equation  (39)  are  also  shown  in  Figure  18.  The 
functions  in  Equation  (39)  to  be  integrated  are  plotted  in  Figure  19  for 
each  wedge  angle. 

The  earlier  analysis  (Reference  3)  of  the  immersing  wedge  was  shown 
to  L.e  in  agreement  with  experimental  data  at  deadrise  angles  of  ap¬ 
proximately  20°.  For  comparison  with  the  present,  more  detailed  analysis, 
the  pressure  distributions  and  loads  from  the  expanding-plate  analogy 
have  been  computed  and  are  plotted  in  Figures  17  and  18  (according  to 
Equation  46  of  Reference  3)  for  20°,  30°,  40°,  and  50°  deadrise. 

In  Equation  76  of  Reference  8,  the  solution  for  total  load  on  an 
18°  wedge  was  given  as* 


P  »  49  8  T  p  V* 


where  T  is  equivalent,  in  the  notation  of  this  report,  to  V  -  (2c/wJtanA. 
Then, 
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P 

c(pV7!2) 


-  49.8  x—  tan  18°  x  2  -  20.7  (two  sides) 

TT 


This  value  for  p  =  18°  is  indicated  in  Figure  18,  along  with  Wagner's  ap¬ 
proximation  for  the  variation  of  load  with  ^eadrise  (Equation  78  of  Ref¬ 
erence  8), 


P 

c(pV*/2) 


=  2  tanA 


/90? 
V  P 


Comparison  of  the  various  load  results  shown  in  Figure  18  indicates 
several  noteworthy  trends.  As  was  to  be  expected  (within  the  accuracy  of 
the  drafting  and  computation  required),  the  loads  obtained  by  the  present 
analysis  are  in  reasonable  agreement  whether  obtained  by  Wagner's  method 
or  by  integration  of  the  bottom  pressures.  The  expanding-plate  analogy 
indicates  higher  pressures  and  loads  than  the  present  analysis,  with  the 
percentage  difference  increasing  with  deadrise  angle  up  to  50°.  The  single 
point  for  load  at  p  =  18°  given  in  Reference  8  lies  far  below  the  curve 
obtained  on  the  basis  of  the  expanding-plate  analogy  (Reference  3).  This 
point  and  Wagner's  approximate  curve  lie  about  8%  below  an  extension  of 
the  curve  of  loads  determined  by  the  present  analysis. 

Thus,  the  variation  of  load  with  deadrise  angle,  as  determined  by 
the  integration  of  the  bottom  pressure,  may  be  approximated  by 

P  /90°  ,\2 

— -  2.16  tan**  (—  -  l) 


DISCUSSION 


SURFACE  SHAPE 

The  theoretical  analysis  for  the  determination  of  the  spray  thick¬ 
ness  was  concluded  with  a  method  for  checking  the  irrotationality  of  the 
flow  field,  using  an  artifice  involving  the  superposition  of  a  known  flow 
(doublet)  and  the  graphical  application  of  Greene’s  theorem  (Equation 
(31)).  In  spite  of  the  great  care  taken  in  drafting  and  computation 
details,  the  spray  thickness  could  not  be  established  as  precisely  as  de¬ 
sired.  The  curve  presented  in  Figure  10  was  selected  as  representative  of 
the  computed  points,  even  though  the  relative  position  of  the  30°  and  40° 
points  indicate  a  possible  flexure  ?n  the  curve. 

After  the  determination  of  the  loads  and  pressures,  another  check 
on  the  spray  thickness  became  apparent.  Wagner’s  analysis  for  the  total 
load  based  upon  a  simultaneous  solution  of  momentum  and  kinetic  energy 
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relationships  holds  only  for  irrotational  potential  flow.  It  was  noted 
previously  that  the  load  was  almost  directly  proportional  to  the  spray 
thickness.  Also,  it  appeared  from  the  actual  calculations  that  the  velo¬ 
city  and  potential  distributions  were  affected  only  slightly  by  the  value 
of  the  spray  thickness.  Thus,  the  agreement  between  the  results  of  the 
earlier  load  analysis  (Equation  (39) »  with  the  integration  of  the  pres¬ 
sures  on  the  wedge  is  a  measure  of  the  accuracy  of  the  spray  thickness. 

Application  of  this  criterion  to  the  results  given  in  Figure  18  in¬ 
dicates  that  the  spray  thicknesses  for  4G°  and  50°  deadrise  are  too  large, 
while  at  20°,  the  estimated  spray  must  have  been  too  thin.  If  the  load 
analyses  are  assumed  to  be  accurate,  then  the  revisions  in  spray  thick¬ 
nesses  required  to  give  consistent  force  resiilts  are  proportional  to  the 
difference  between  the  computed  load  points^A  revised  set  of  soray  thick¬ 
ness  points  was  determined  in  this  manner  and  plotted  in  Figure  10. 
These  revised  values  of  the  spray  thickness  are  very  well  fitted  by 


(60) 


It  should  be  noted  that  the  spray  thickness  is  of  a  small  order  of 
magnitude.  At  20°  deadrise,  the  total  thickness  is  only  slightly  more 
than  1 %  of  the  we" ted  half-width,  while  atthe  highest  angle  investigated, 
the  spray  thickness  ratio  is  only  8%.  (This  factor  will  make  experimental 
verification  of  the  surface  o'.ape  very  difficult  except  in  a  qualitative 
way.)  In  view  of  the  amount  of  graphical  construction  involved  in  the 
present  analysis,  the  spread  between  the  spray  thickness  curves  in  Figure 
10  is  not  unreasonable. 

If  extension  of  the  computations  to  higher  deadrise  angles  becomes 
desirable,  it  will  probably  be  sufficient  to  by-pass  the  doublet  check  or. 
irrotationality  and  go  directly  to  the  determination  of  loads  and  pres¬ 
sures  for  surface  curves  based  on  the  extrapolated  values  of  5/c  from 
Equation  (60).  Since  the  computation  of  loads  and  pressures  is  not  too 
tedious  for  the  estimated  spray  thickness  and  surface  shape,  and  since 
the  approximate  correction  to  the  spray  may  be  estimated  using  the  ap¬ 
proximate  proportionality  between  the  load  and  spray  thickness,  it  would 
be  economical  to  use  the  lead  integrals  as  the  check  on  the  irrotation¬ 
ality  of  the  flow  field. 

The  emphasis  in  this  discussion  has  been  upon  the  determination  of 
the  relative  spray  thickness.  This  is  so  because  the  shape  of  the  surface 
may  be  very  closely  approximated  bv  the  graphical  construction  method  for 
similiarity  and  continuity  once  the  proper  spray  thickness  is  knotyn. 


VELOCITY  AND  POTENTIAL  DISTRIBUTIONS  ALONG  THE  HEDGE 


It  has  been  noted  that  the  velocity  along  the  wedge  as  determined 
from  this  analysis  is  in  partial  agreement  with  that  indicated  by  the  ex¬ 
panding-plate  analysis.  As  discussed  in  Reference  3,  the  plate  analogy 
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breaks  down  near  the  edge  with  the  indication  of  infinite  velocity.  That 
would  be  a  physical  impossibility,  and  in  actuality,  a  spray  root  is 
formed.  The  velocity  increases  less  rapidly  than  indicated  by  the  expand¬ 
ing-plate  analysis,  as  indicated  in  Figure  16,  and  then  levels  off  to  a 
constant  value  in  the  spray. 

By  virtue  of  the  symmetry,  at  the  centerline  of  a  flat  expanding 
plate  the  velocity  along  the  plate  must  be  zero.  However,  on  the  actual 
wedge,  there  is  a  component  along  the  wedge  toward  the  center  so  that  the 
resultant  velocity  at  the  apex  is  vertical  (not  normal  to  the  bottom  of 
the  wedge). 


This  is  the  difference  which  appears  at  x/c  -  0  in  Figure  16.  The 
close  agreement  between  the  velocity  curves  in  the  range  0.4  <  x/c  <  0.95 
is  quite  remarkable,  expecially  for  the  30°  deadrise  wedge.  The  increase 
with  increasing  deadrise  in  the  difference  between  the  present  curves 
and  those  from  the  expanding-plate  analogy  appears  to  be  due  primarily  to 
the  condition  at  the  step  apex  discussed  above. 


With  respect  to  the  potential  on  the  wedge  as  given  by  the  present 
analysis  compared  to  that  from  the  expanding-plate  analogy,  a  more  in¬ 
volved  reasoning  is  required.  It  will  be  noted  in  Figure  14  that  the  pre¬ 
vious  and  present  analyses  give  the  same  shape  of  potential  distribution, 
but  that,  the  expanding-plate  analogy  indicates  higher  negative  values 
(almost  a  constant  difference  along  s/c  for  any  particular  deadrise 


angle).  I'uuet  the  assumptions  of  Deference  3  (expanding— plate  analogy), 


the  potential  is  based  at  zero  on  the  level  of  the  plate.  This  becomes  an 


untenable  position  at  any  appreciable  deadrise  angle,  since  the  level 


water  surface  lies  above  the  wedge  apex,  and  the  plate  cannot  be  construed 


as  being  at  several  levels  simultaneously.  The  difference  in  potential  of 
the  field  for  a  reference -level  change  of  yQ  (=(2c/?0 tan^fl  ]  is  VyQ.  In 
terms  of  the  scale  on  Figure  14,  the  change  in  7  is 
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This  is  almost  exactly  the  difference  at  the  lower  deadrise  angle&;between 
the  expanding-plate  potential  and  the  distribution  given  by  the  present, 
more  exact  analysis. 


PRESSURES  AND  LOADS 

Because  of  -the  -independent  check  of  the  spray  thickness  by  the  com¬ 
parison  of  total  load  and  integrated  pressure  distributions,  it  may  be 
reasoned  that  the  pressures  are  reasonably  accurate.  This  could  not  be 
assumed  to  be  the  case  if  such  a  check  had  not  been  obtained,  because  of 
the  number  of  tpsrms  and  graphical  constructions  involved  in  the  steps 
leading  to  the  determination  of  the  pressure. 

The  failure  of  the  expanding-plate  analogy  in  predicting  the  correct 
wedge  pressures  may  be  partly  due  to  the  difference  in  potential  which 
was  discussed  above.  It  may  be  shown  (by  reference  to  Equation  (51))  that 
a  change  in  potential  (&£  =  Fy0)  shows  up  in  the  pressure  ratio,  P/(pV2/ 2), 
as  an  additive  term  (of  value  -2).  Thus,  a  changed  potential  reference  in 
the  expanding-plate  analogy  will  reduce  the  pressures  shown  in  Figure  17 
by  2.  This  does  not  bring  all  of  the  curves  into  alignment,  but  a  trend 
in  that  direction  is  obviously  correct. 

It  may  be  seen  in  Figure  17  that  the  trends  of  tta  pressure  differ¬ 
ences  at  the  peak  (in  the  spray  root)  and  at  the  keel  are  opposite.  Al¬ 
though  the  expanding-plate  maximum  pressures  are  higher  than  those  ob¬ 
tained  from  the  present  detailed  analysis,  the  percentage  excess  decreases 
with  decreasing  deadrise.  On  the  other  hand,  the  keel  pressure  given  by 
the  earlier  analysis  is  higher  at  low  deadrise  angles  (approximately  by 
the  amount  of  2,  as  mentioned  above),  but  the  difference  decreases  with 
increasing  angle. 

It  is  interesting  to  note  that  a  rather  uniform  pressure  distribu¬ 
tion  is  obtained  at  40°  deadrise  angle.  At  lower  angles,  the  peak  pres¬ 
sures  greatly  exceed  the  average  value  across  the  section,  while  at  higher 
angles,  the  lack  of  a  well-defined  stagnation  area  in  the  spray-root 
region  lets  the  pressure  drop  below  the  value  at  the  keel.  This  trend  had 
been  shown  before  in  the  expanding-plate  analysis,  but  could  not  be  ac¬ 
cepted  upon  that  basis  alone  because  of  the  wide  divergence  from  the  con¬ 
ditions  assumed  for  that  analysis.  However,  the  removal  of  the  restriction 
to  low  angles  made  possible  by  the  present  analysis  has  resulted  iu  a 
similar  variation  of  pressure  with  deadrise  angle,  so  that  the  trend  may 
be  considered  to  be  established, 


CONCLUSIONS 


1.  A  method  has  been  presented  for  the  theoretical  determination  of 
the  free  surface  shape  arising  from  the  symmetrical  penetration  of  the 
surface  of  a  perfect  fluid  by  an  infinite  wedge-shaped  body  at  constant 
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velocity  (normal  to  the  initial  undisturbed  surface) 

2,  The  actual  surface  shapes  have  been  constructed  for  20°,  30°, 
40°,  and  50°  deadrise  angles. 

3.  Since  no  restriction  upon  the  wedge  angle  or  slope  of  the  free 
surface  was  necessary  to  this  solution,  the  method  may  be  used  for  even 
greater  deadrise  angles  when  desirable  . 

4  In  the  process  for  constructing  these  surface  shapes,  the  spray 
thickness  ratio  appeared  as  the  most  critical  dimension.  The  spray  is 
uniformly  tapered  from  zero  thickness  at  the  tip  to  a  maximum  in  the 
spray  root  where  it  is  related  to  the  wetted  half-width  by  the  empirical 
relationship  $/c  =  (>3/ir)2, 

5.  Distributions  of  the  velocity,  potential,  and  pressure  along  the 
wedge  surface  were  evaluated  by  the  extension  of  the  solution  for  the 
free  surface  shape  These  distributions  are  similar  to  those  given  by  the 
earlier  expanding-plate  analysis  (Reference  3)  except  that  an  accurate 
description  of  the  flow  is  now  given  in  the  spray-root  region, 

6  The  total  loads,  obtained  by  the  integration  of  the  pressure 
over  the  wetted  width  of  the  wedge,  are  given  by  the  empirical  relation¬ 
ship 


P 

c(pV*I2) 


-  2  16  tan 


I 


which  is  8%  higher  than  Wagner's  approximation  of  the  same  form. 
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APPENDIX  I 

SAMPLE  CONSTRUCTION  OF  THE  FREE  SURFACE  SHAPE 

This  appendix  presents  the  details  of  the  computations  and  drafting 
methods  employed  in  the  construction  of  the  free  surface  shape  that  meets 
the  conditions  of  continuity,  similarity,  and  irrotationality.  To  obtain 
greater  clarity,  not  all  of  the  construction  lines  are  drawn  and  some  of 
the  tables  are  abbreviated. 


FIRST  ESTIMATE  OF  THE  SURFACE  SHAPE 

A s  noted  in  the  body  of  the  report,  the  free  surface  is  built  up  of 
spray  and  wave.  The  first  estimate  of  the  wave  rise  is  given  by  Equation 
(3).  For  this  particular  case,  c  is  taken  as  10  inches  for  A  =  40°  and 


FIGURE  1-1 


The  spray  tip,  o',  is  located  by  measuring  along  the  wave  and  wedge 
from  x  -  15  a  distance  s  -  x  =  15.  (At  the  distance  x  ~15{  the  wave 
slope  is  very  small  SO  that  i  -  x  is  a  very  close  approximation. }  A  splay 
thickness  r-tio,  S/c  =  0.055,  was  assumed. 

The  addition  of  spray  requires  that  the  wave  rise  be  reduced  to 
satisfy  continuity  of  the  flow.  In  order  to  keep  a  simple  analytic  form 
for  the  wave  at  a  distance,  the  value  of  v0  is  reduced  proportionately: 

Vj  =  «n0 


36-- 


and  the  spray  is  blended  into  the  reduced  wave  to  form  the  spray-root  re 
gion. 


The  spray  tip  is  located  more  exactly  by  following  the  curved  sur¬ 
face,  and  area  B  plus  C  is  made  to  equal  area  A  by  varying  the  wave  rise 
ij  / 1)  and  to  a  small  extent,  by  the  spray-root  curve.  The  area  C  beyond 
x  15  on  the  wave-rise  curve  may  be  obtained  by  the  integral  of  Equation 
(3): 


7]  dx  -  a  j 7)0  dx  (see  Equation  (4)). 

^  w 


APPLICATION  OF  SIMILARITY  CONDITION 


The  estimated  surface  shape,  as  drawn  above,  is  of  the  correct 
length  and  encloses  the  proper  fluid  area.  It  is  then  necessary  to  apply 
the  test  of  Equation  (17),  using  the  assumed  shape  to  determine  the  velo¬ 
city  components  (Figure  3).  The  actual  calculations  for  this  particular 
case  are  given  in  Table  I. 


If  the  first  approximation  to  the  surface  shape  is  in  error,  the 
vaJues  of  n  and  f  resulting  from  this  graphical  integration  wi.l  differ 
from  the  assumed  values.  Then,  the  process  is  repeated,  using  the  newly 
calculated  coordinates  of  the  free  surface  until  a  close  check  is  obtained 
in  Table  1  between  the  starting  and  final  v  ana  f. 
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APPL i CAT  I  ON  OF  CONDITION  FOR  I RROTAT I ONAL  FLOW 

The  superposition  of  the  potential  flow  due.  to  a  doublet  upon  the 
flow  pattern  determined  by  the  free  surface  shape  obtained  above  must 
satisfy  Equation  (31)  for  irrotational  flow.  The  quantities  t4>  and 
t  O^/3/i )  for  the  trial  surface  shape  are  obtained  by  Equations  (20)  and 
(15)  as  shown  in  Figure  I-l  (resolving  tU  parallel  and  normal  to  the  sur¬ 
face).  The  valuts  of  and  for  the  doublet  are  given  by  Equations 

(25)  and  (28)  and  are  computed  in  Table  II.  The  values  of  t<£  and  t(3?fr/3n) 
are  computed  in  Table  III. 

From  Figure  I-l,  the  left-hand  side  of  Equation  (31)  (multiplied  by 
t)  may  be  evaluated.  For  the  particular  case  of  $/c  =  0.055  for  /3=40° 
(c  =  10  inches),  the  value  of  t  ^1  at  the  doublet  (with  respect  to  the 
doublet)  is  s'^ 


sd  ~  Jo  C0SI 3  =  ‘6.18 


Also.  tV  cos/3  =  y0  cosfi  =  4.09  and  p  -  2.86  radians.  Thus, 


B*  „  -(-6.18)  2  (4.09K2.86) 

—  *  2V/i  cos/3  =  -  +  - 

**  d  *  * 


* 42.8 

t 


The  right-hand  side  of  Equation  (31)  is  evaluated  by  computing  the 
values  of  [t<2>(  3<^/3n )  -  <t>^  t(3tf/3n)3  and  integrating  along  the  surface. 
This  is  done  in  Table  III,  making  use  of  the  values  of  tfy  and  3^/3n  from 
Table  II-  In  this  particular  case,  an  approximately  correct  spray  thick¬ 
ness  was  chosen  so  that  value  of  the  integral 


In  general,  however,  several  spray  thicknesses  would  be  assumed  and 
the  complete  process  'carried  out  to  this  stage.  Then  plots  of  the  left- 
and  right-hand  sides  of  Equation  (31)  vs.  i/c  would  determine  by  their 
crossing  the  correct  value  of  the  spray  thickness  ratio. 
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TABULAR  INTEGRATION  FOR  SURFACE  SHAPE  (EQUATION  (17)) 
(.0  -  40°,  c  =  10  inches,  5  5  0.55  inches) 


mm 

tv 

tv 

Jtv 

tu 

(U 

tu 

5 

tv 

-As 

V 

tu 

-2 

3As 

l~As 

f -*o 

& 

s 

s 

NEE 

START 

s 

r 

3 

NEW 

15.0 

0  37 

1.43 

0.0064 

0.018 

0.0247 

0.37 

13.0 

0.56 

2.13 

0.018 

0,0427 

0.56 

0 

0.09 

0.0005 

0.001 

0 

0 

12.0 

0  72 

3.21 

0.0224 

0.014 

0.0607 

0.73 

0.01 

0.26 

0.0018 

0.0015 

0.001 

0.012 

11.5 

0.86 

4.12 

0.0312 

0.0185 

0,0747 

0.86 

0.03 

0.50 

0.0038 

0.0030 

0.0025 

0.029 

11  0 

1  02 

5  45 

0.0450 

9*9212 

0.0932 

1.025 

0.06 

1.00 

0  0083 

0.0644 

0.0055 

0  060 

1C. 6 

1  21 

7  0! 

0.0624 

-J* 

0.1144 

1.21 

0.11 

1.88 

0.0167 

* 

0.0099 

0.1  OS 

0.0386 

0*0116 

1C. 2 

1.48 

9.87 

0.0949 

0.0432 

0.1450 

1.48 

0.22 

4.64 

0.  0446 

0.0276 

0.021S 

0.219 

9.8 

i.ts 

11.61 

0.1210 

0.9496 

0.1882 

1.85 

0  48 

9.28 

0.0967 

0.0508 

0.0491 

0.481 

HI 

2  24 

a.  o: 

0.1252 

0.0502 

0.2378 

2.24 

0.94 

13.60' 

8.1540 

0.0692 

0.0999 

0.939 

9.0 

2  S9 

>0.16 

0.125S 

0.064 

0.2880 

2.59 

1.52 

15.41 

0.1903 

0.1045 

0.1691 

1.521 

8.S 

2  99 

9  42 

0.1303 

0.C663 

0.3520 

2.99 

2.33 

16.40 

0.2272 

0.12! 

0.2736 

2.325 

8.0 

3.36 

9.16 

0  1431 

0.161 

0.4203 

3.36 

3.16 

16.68 

0.2605 

0.3025 

0.3946 

3.165 

7.0 

4  07 

8  87 

0.1810 

0.5813 

4.07 

4.88 

16.98 

0.346 

0.6971 

4.88 

TABLE  II 

COMPUTATION  OF  $d  AND  2<j>d/2n 
=  40°,  c  =  10  inches,  5  =  0.55  inches) 
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9.67 
9.87 
tO.  07 
10.27 
10.47 
1  0.67 
1  0.87 
11.0 
11. S 
12.0 
13.0 


1 

.16 

1 

.82 

0 

.58 

0 

.39 

0 

.20 

0 

0 

.20 

.39 

1 

.59 

71 

169.0 
165.  S 
158. 0 

149.5 
143.0 

132.5 
121.0 
111.0 

102. 5 
96. 5 
90.0 

96.5 

99.5 

192. 5 

105.5 

112.5 
116.0 
121.0 


166.0 

159.0 

145.0 

131.0 

121.5 

109.5 

100.0 

9S.0 

92.0 

90.0 

90.0 

90.5 
91.0 
92.0 
93.0 

96.5 

99.5 
104.5 


6d 

^ d 

an 

e.o 

17.5 

0.063 

0.  002 

11.5 

24.5 

0.121 

0.010 

18.5  i 

38.5 

0.294 

0.  060 

24.5 

54.5 

0.S27 

0.167 

28.5 

67.5 

1 

0.724 

0.260 

30.5 

88.0 

1.039 

0.324 

29.5 

109. S 

1.455 

0.210 

23.5 

130.5 

2.004 

0.757 

15.0 

150.0 

2.876 

0.476 

7.0 

165.5 

5.318 

0.249 

0.3274® 

-4.661 

-2.239 

-1.381 

-1.091 

-0.S28 

-0.317 

-0.144 


0 

-0.671 
-0.137 
-0.148 
-0.180 
-0.156 
-0.  099 
-0.054 


FIGURE  2 
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GENERAL  SHAPE  OF  THE  FREE  SURFACE  OF  A  FLUID 
DUE  TO  PENETRATION  BY  A  WEDGE 
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figure  4 


LOCATION  OF  DOUBLETS  FOR  I  RROTATiONALITY  CHECK 
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+  $  df  =  $  cos/?  ds 


FIGURE  6 

SIGN  CONVENTIONS  FOR  SURFACE  ELEMENTS  AND  VELOCITIES 


FIGURE  7 

POLAR  COORDINATE  SYSTEM  FOP  SOLUTION 
OF  CAUCHY  INTEGRAL  EQUATION 
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COMPARISON  OF  FIXED  AND  RELATIVE  POSITION  POINTS 

IN  THE  FLUID 


FIGURE  I! 


SUPERPOSITION  OF  SEVERAL  SPRAY-  ROOT  SOLUTIONS 


VELOCITY  DISTRIBUTION  ON  WEDGE  AND  FREE  SURFACE 


VELOCITY  DIRECTION  AT  WEDGE  AND  FREE  SURFACE 


FIGURE  17 


FIGURE  18 
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